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ABSTRACT 



The 09V star HD 46202, which is a member of the young open cluster NGC 2244, was observed by the CoRoT satellite in 
October/November 2008 during a short run of 34 days. From the very high-precision light curve, we clearly detect /? Cep-like pulsation 
frequencies with amplitudes of ~0. 1 mmag and below. A comparison with stellar models was performed using a.x^ as a measure for 
the goodness-of-fit between the observed and theoretically computed frequencies. The physical parameters of our best-fitting models 
are compatible with the ones deduced spectroscopically. A core overshooting parameter o-o,, = 0. 10 ± 0.05 pressure scale height is 
required. None of the observed frequencies are theoretically excited with the input physics used in our study. More theoretical work is 
thus needed to overcome this shortcoming in how we understand the excitation mechanism of pulsation modes in such a massive star 
A similar excitation problem has also been encountered for certain pulsation modes in /3 Cep stars recently modelled asteroseismically. 

Key words. Asteroseismology - Stars; oscillations - Stars: variables: early-type - Stars: individual: HD 46202 



1. Introduction 

Asteroseismology is the domain of astrophysics that probes 
the internal structure of stars by using their stellar oscillations. 
Different classes of pulsating stars are known across the H-R 
diagram, including pulsating main-sequence B-type stars (with 
masses between 3 and 20 Mq), such as the so-called /3 Cep 
stars and the slowly pulsating B stars (e.g., Aerts et al. 1 20101 ). 
Currently, the part of the H-R diagram corresponding to even 
hotter and more massive main-sequence stars has been poorly 
explored for asteroseismic purposes, both theoretically and ob- 
servationally. 

Due to lack of firm observational detections, the theory of 
stellar oscillations has hardly been applied to O star models. 
On the one hand, one expects, for late O-type stars such as 
HD 46202 (09V) studied in this paper, that oscillation modes 
excited by the k mechanism occur with periods of several hours, 
just as for the well-established class of /3 Cep stars. The theo- 
retical treatment for the lower mass main-sequence O stars is 
the same as the one for the early-type B stars because stel- 
lar winds are not expected to play an important role in the os- 
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dilatory behaviour. Radiation-driven winds have thus been ig- 
nored when computing instability strips (e.g., Pamvatnvkh l 19991 
Miglio et al. 2007 a b). On the other hand, theory predicts the oc- 
currence of so-called strange mode oscillations in stars close to 
the Humphreys-Davidson limit. These are due to a convection 
zone in the outer stellar layers where a density inversion occurs 
and the radiation pressure dominates (e.g., Glatzel & Kiriakidis 
11993, Papaloizou et al. ,1997. Saio et al. 119981) . Review papers 
on the topics are available by Glatzel (120091 ) and by Saio (2009|. 
Observations have not established that there are strange modes 
in O stars. 

On the observational side, there are a few O-type stars for 
which spectroscopic and photometric variability linked to stellar 
pulsation has been searched for and found from ground-based 
data. 

Convincing evidence of stellar pulsation has been detected in 
( Oph (Kambe et aL [T997] ), HD 93521 (Howarth & Reid 1993] 
Rauw et al.|2p08 ), ^ Persei, and A Cephei (de Jong et al. 199^. 
In fact, it is difficult to discover pulsation modes in O-type tar- 
gets because the observed variability that might come from pul- 
sations has very low amplitude (e.g., Henrichs 1 1 9991 for a sum- 
mary). 

In this context clearer detections can be expected from long 
uninterrupted photometric time series at micro-magnitude preci- 
sion obtained from space satellites such as CoRoT (Convection, 
Rotation and planetary Transits, Auvergne et al. 120091 ) . As part 
of the asteroseismology programme, several O-type stars were 
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observed by CoRoT during the second short run SRa02 point- 
ing towards the anticentre of the Galaxy, e.g., the O-star bina- 
ries HD 46149 (Degroote et al. 2010) and HD47129 (Mahy et 
al. [2010). In this paper, we present a study of the single 09V 
star HD 46202, which is a member of the young open cluster 
NGC 2244 in the Rosette Nebula. 



2. Frequency analysis on CoRoT photometry 

The CoRoT data of HD 46202 were taken during 34 days be- 
tween 2008 October 8 and 2008 November 12 with an average 
time sampling of 32 seconds. This involves a frequency reso- 
lution of ~0.03d"'. After removing all flagged measurements, 
such as the ones suffering from hot pixels during the passage 
through the South Atlantic Anomaly, the light curve consists 
of 81 781 datapoints. As in most CoRoT light curves, a down- 
ward trend considered to be of instrumental origin is visible. 
Moreover, the light curve presents three jumps during the last 
five days of observation. We divided the light curve into four 
parts and, for each of these selected regions, we removed a trend 
using a linear approximation. To visualise the variability, only 
part of the detrended light curve is shown in Fig.[T] We can 
clearly see a variation with a period of the order of 2 days. 
Moreover, variations on a shorter time scale are also observed. 

The Scargle periodogram (Scargle 119821) of the detrended 
light curve is presented in Fig. [J] Clear peaks are detected be- 
tween 0.5 and 5 d ' . The peak at 13.97 d"' is not intrinsic to the 
star but corresponds to the orbital frequency of the CoRoT satel- 
lite. To determine the whole frequency spectrum of the star we 
used the methodology described in Degroote et al. (2009a), who 
analysed numerous B-type pulsators using light curves from the 
exoplanet data of the CoRoT satellite. This methodology uses 
the usual method of successive prewhitening, along with a statis- 
tical stop criterion based on the one first introduced by Scargle 
( 119821 ). The statistical stop criterion defined in Degroote et al. 
( |2009al ) was adopted instead of a signal-to-noise (S/N) criterion 
(Breger et al. 119931 ). The reason is that, although a signal origi- 
nating from nonsinusoidal processes (e.g., red noise) is also con- 
sidered to be significant, and thus prewhitened using a sum-of- 
sines model, isolated peaks at lower amplitudes can be recovered 
outside of these regions. 

The frequencies with the highest amplitudes, which turn out 
to be typical /JCep-like pulsation modes, are listed in Table[T] 
In addition, all the frequencies with an amplitude lower than the 
one corresponding to the satellite's orbital frequency are given in 
the Appendix. In total 67 frequencies are found to be significant 
by our adopted stop criterion. For comparison purposes, we also 
indicate the signal-to-noise (S/N) ratio, which is computed using 
a 6d ' box-interval centred on the considered frequency. The 
error estimate on the frequency is computed as in Degroote et al. 
( |200 9a^, i.e., using the formula of Montgomery & O'Donoghue 
( 1199 9 ). but also taking correlation effects into account according 
to the formalism by Schwarzenberg-Czerny (l2003b . 

Recent studies of CoRoT light curves of B-type stars high- 
lighted variations in the amplitude of frequencies and in the fre- 
quencies themselves as a function of time (e.g., in the jSCep 
star HD 180642, Degroote et al. 2009bi in the BO.SIVe star 
HD 49330, Huat et al. 2009). In order to check the stabiHty 
of the frequencies of HD 46202 listed in Table[Tl we applied 
a short time Fourier transform. Different window widths were 
used, leading to the conclusion that the frequencies are stable in 
time. There is no convincing proof of changing amplitudes ei- 
ther. The detected low amplitude variations in time are typical of 
beating effects for a multiperiodic pulsator (Fig. [3]). 
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Fig. 1. Part of the detrended CoRoT light curve, converted to 
magnitude (mag - -2.5 log[o(flux), once the average has been 
substracted). The first datapoint corresponds to 2454759 HJD. 
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Fig. 2. Scargle periodogram of the detrended CoRoT light curve, 
converted to magnitude. 



We searched for rotational frequency splittings with the au- 
tomatic methods described in Degroote et al. (2009a), but we 
could not unambiguously or confidently find such multiplets. For 
the highest amplitude modes, the frequencies /4, /s, and (see 
Table[TJ might belong to the same multiplet and the frequencies 
/g, /lo, and /ii might belong to another one. However, recent 
work on the fi Cep star 12 Lacertae taught us to be cautious in 
attributing frequencies to multiplets without empirical identifi- 
cation of the mode degree. For this latter target, Dziembowski & 
Jerzykiewicz ( 1999) assumed that three of the strongest modes, 
almost equidistant, belong to the same multiplet in their seismic 
modelling. Later, a photometric mode identification by Handler 
et al. (120061 ) ruled out this assumption by proving that the three 
modes are actually associated with three different values of the 
degree I. It is obvious that reliable empirical mode identification 
from multicolour photometry, and/or spectroscopy, is indispens- 
able for definite conclusions. 

CoRoT observations also revealed solar-like oscillations in 
two massive stars: the Cep star HD 180642 (Belkacem et al. 
2009) and the O-type binary HD 46149 (Degroote et al. 201^. 
To search for the signature of such oscillations in HD 46202, 
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Table 1. The highest-amplitude frequencies derived from the de- 
trended CoRoT light curve after subsequent prewhitening, their 
corresponding amplitude, and signal-to-noise ratio. 



ID 






Amplitude (mmag) 


S/N ratio 




0.510 


5.901 


0.1079 


9.3 


fi 


4.856 


56.204 


0.0934 


11.7 


h 


4.691 


54.295 


0.0798 


10.7 


/4 


2.290 


26.500 


0.0708 


8.5 


h 


2.643 


30.586 


0.0645 


7.7 


h 


3.004 


34.769 


0.0587 


7.2 


fi 


2.195 


25.404 


0.0460 


5.7 


h 


1.462 


16.919 


0.0426 


5.6 


h 


2.222 


25.722 


0.0339 


4.8 


/lO 


0.804 


9.305 


0.0314 


4.5 


/ii 


2.127 


24.613 


0.0307 


4.3 


f'n 


13.973 


161.728 


0.0298 


10.6 
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Fig. 3. The short time Fourier transform of the detrended CoRoT 
light curve, computed with a window width of 10 days. The 
colour scale is in mmag units. 
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Fig. 4. Temporal variance spectrum (TVS) computed from the 
cross-correlated profiles computed by taking all elements into 
account, except H. The line corresponds to the 99% signifi- 
cance level for the variability evaluated following the approach 
of Fullerton et al. ( fT996t . 



First, variability was searched in individual line profiles. Then, 
to increase the S/N of the spectra and improve the detection 
of the line-profile variations, we used cross-correlated profiles 
computed using the least squares deconvolution (LSD) method 
(Donati et al. 11997 ) and taking all elements into account except 
H. We also computed cross-correlated profiles associated with 
only one element for all possible elements, such as helium, sil- 
icon, nitrogen and oxygen. To perform our frequency analysis, 
we used the software package FAMIASQ (Frequency Analysis 
and Mode Identification for Asteroseismology, Zima 2008). The 
significance of variability in stellar spectra can be quantified by 
a temporal variance spectrum (TVS, Fullerton et al. 1996). We 
computed TVS diagrams for individual line profiles and cross- 
correlated profiles (Rauw et al. l2001l l. One of them is illustrated 
in Fig.|4] Even with such state-of-the-art techniques, no spectro- 
scopic variability could be detected. 



we performed an autocorrelation of the time series and peri- 
odogram, both for the original detrended light curve and for the 
light curve prewhitened with the highest-amplitude frequencies. 
We could not find any structure typical of solar-like oscillations 
with amplitudes above 0.009 mmag. 

3. Search for variability in spectra 

From the photometric amplitudes of the pulsation modes being 
0.1 mmag or less, it is clear that these modes cannot be de- 
tectable in spectra of moderate S/N. Nevertheless, the study of a 
time series of spectral lines remains relevant when searching for 
other kinds of variability. First, we can check for the presence of 
a variable stellar wind. Also, modes of very high-degree t might 
be visible in the spectral lines (e.g., as in the ySCep star <y' Sco, 
Telting & Schrijvers 1998), although they are not observed in 
the CoRoT light curve. Additionally, some spectral lines associ- 
ated to only certain elements might be variable with a periodicity 
linked to the stellar rotation and magnetic field. 

Consequently, we gathered 44 high-resolution spectra (R = 
50000, S/N of -60) with the CORALIE echelle spectrograph 
attached to the 1.2m Leonard Euler telescope (La Silla, Chile). 



4. Fundamental parameters 

In the next section, we compute stellar models for our target, 
which match the pulsation frequencies discovered in the CoRoT 
photometry. The fundamental parameters and chemical compo- 
sition of HD 46202 are also used as additional constraints. To 
deduce the latter, we used a high-quality spectrum with ex- 
tended wavelength coverage, taken with the ESO/MPG 2.2m 
telescope with FEROS (Kaufer et al. 1999). Six individual spec- 
tra, observed between 2008 April 16-26 of 720 s exposure time 
each were extracted from the ESO archive and reduced with the 
MIDAS pipeline. Coaddition of these data gave a spectrum with 
S IN of more than 500 per resolution element in the optical blue 
region. 

Our quantitative spectral analysis follows the methodology 
described by Nieva & Przybilla (l2007l 12008]) . In brief, non- 
LTE line-formation calculations were performed with the codes 
Detail and Surface (Giddings TMT: Butler & Giddings [19851 
both updated by K. Butler) on the basis of a prescribed LTE 
atmospheric structure (Atlas9, Kurucz I1993I I and accounting 

' FAMIAS has been developed in the framework of the FP6 European 
Coordination Action HELAS ,http://www.helas-eu.org/ 1 
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Fig. 5. Examples of line fits (red) for H, He, 
and C to the high-quality observed spectrum 
(black). The hydrogen and helium lines are 
reproduced with a single synthetic spectrum, 
while for carbon a slightly different abundance 
value is adopted per line, giving rise to a 1-cr 
uncertainty of 0. lOdex due to line-to-line scat- 
ter. 



Table 2. Atmospheric parameters, photometric properties and fundamental parameters of HD 46202 as derived from the spectro- 
scopic analysis. 



reff(K) 34100+600 V 8.18±0.02 M/Mq 19.6±2.7 

logg(cgs) 4.17±0.07 B-V 0.17+0.02 R/R^ 6.5+_°-^^ 

^(kms-') 6±2 E{B-Y) 0.49±0.03 logL/Lo 4.71 ±0.10 

vsini(kms-') 25±7 Ry 2.7±0.1 X (by mass) 0.715±0.033 

^(kms-i) 15±7 Mv -3.73±0.25 y(bymass) 0.271±0.033 

M(Moyr-i) <10-' B.C. -3.30+0.05 Z(bymass) 0.014+0.001 

rfspec(kpc) 1.3+0.2 Mboi -7.03+0.26 t^voi (Myr) 2+1.5 



for the most recent model atoms of Przybilla et al. ( 120081) . The 
method is equivalent to full non-LTE model atmosphere analyses 
for late O-type stars with weak winds (Nieva & Przvbilla l2007l l. 

The parameters were determined via simultaneous fitting 
of all observed H lines and ionisation equilibria of He i/ii and 
C m/iv. Our results from the spectral fitting are summarised in 
Table 12] for the effective temperature Teff, surface gravity log^, 
microturbulence ^, and the projected rotational v sin / and macro- 
turbulent velocity (. Examples of our fits to the observed spec- 
trum are shown in Fig. |5] for these final parameters. Additional 
models were calculated with the non-LTE code Fastwind (Puis 
et al. I2005I I. which accounts for non-LTE metal line blanketing, 
mass loss, and spherical geometry, for an independent test of the 
parameter determination. Similar atmospheric parameters were 
derived from observed H and Hei/ii line spectraQ, confirming 
the match of our hydrostatic approach also with hydrodynam- 
ical models for late O-type main sequence stars with weak stel- 
lar winds. The Fastwind calculations also allowed the mass-loss 
rate M from Ha line-profile fits to be constrained. The upper 
limit found can be considered as insignificant in our comparison 
with stellar evolutionary models. 

The helium, carbon, nitrogen, and magnesium abundances 
derived with Detail/Surface coincide with the cosmic abun- 
dance standard (CAS) from early B-type stars in the solar neigh- 
bourhood as proposed by Przybilla et al. (2008). From an exten- 
sion of that work (Nieva & Przybilla, in preparation) we found 
evidence that abundance analyses of O, Si, and Fe for late O-type 
stars require extensions of the model atoms, which is beyond the 
scope of the present work. We therefore adopt CAS values for 



^ A comprehensive analysis of the metal line spectrum is not feasible 
with Fastwind at present. 



these elements, implying mass fractions for hydrogen X, helium 
y, and the metals Z in HD 46202 as indicated in Table |2] impli- 
citly assuming that the present surface abundances are represen- 
tative for the pristine values. 

To verify the results from our spectroscopic analysis, we 
compared our non-LTE model flux with the observed spectral 
energy distribution of HD 46202 in Fig. |6] For this UV spec- 
trophotometry was extracted from the lUE archive, and Johnson 
magnitudes were adopted from Morel & Magnenat (1978). The 
latter were transformed into absolute fluxes using zeropoints 
from Bessel et al. (1998). The observed flux was dereddened ac- 
cording to the reddening law described by Cardelli et al. ( 19891), 
with both the colour excess E{B - V) and the ratio of total-to- 
selective extinction Ry requiring determination. The data were 
normalised in V. An excellent match of theory with observation 
is obtained for the values given in Table |2] The strong Lye line 
mostly has interstellar origins. 

This allowed the dereddened apparent visual magnitude to 
be determined and the spectroscopic distance t/jpec to be derived 
following Ramspeck et al. (2001). The required stellar mass M 
was obtained by comparison of the star's position in the log T^s- 
log g diagram with evolutionary tracks computed as described in 
Sect. 5 (also the evolutionary age Tevoi was estimated this way). 
This in turn facilitated the absolute visual magnitude My to be 
determined, along with the bolometric magnitude Mboi, the stel- 
lar luminosity L, and the stellar radius R, by using the bolo- 
metric correction B.C. from our model. Thus, a fully consistent 
spectroscopic solution for the fundamental stellar parameters of 
HD 46202 was established, as summarised in Table |2] 

In consequence, HD 46202 is an unevolved star close to the 
zero-age main sequence, which shows no signature of mixing 
with CN-processed material (Przybilla et al. 120101 ). The derived 
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Fig. 6. Fit of the Detail non-LTE model flux for our finafly 
adopted parameters from Table|2](red line) to the observed spec- 
tral energy distribution of HD 46202 (black line, squares) from 
the UV to the near-IR. The observations have been dereddened 
and normalised in V. 
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log Teff 



4.48 



Fig. 7. The 1-cr and 3-cr error boxes represent the position of 
HD 46202, deduced from spectroscopy, in the log Teff-logg di- 
agram. The positions of all the asteroseismic models having a 
<2 are also shown, for X = 0.682 (in red), X = 0.715 (in 
green), and X - 0.748 (in blue). We refer to the text for more 
explanations. 



colour excess, Ry, distance, and age are in excellent agreement 
with the NGC 2244 cluster properties determined by Hensberge 
et al. (I2000l i. 



5. Comparison with stellar models 

We checked that state-of-the-art stellar models with standard 
physics can account for the frequency spectrum observed by 
CoRoT. When X and Z are fixed, we attempted to asteroseis- 
mically constrain the model parameters of HD 46202, which are 
the mass, the central hydrogen abundance (or age), and the core 
convective overshooting parameter. 

5. 1 . Numerical tools and model input physics 

Stellar models for non-rotating stars were computed with the 
evolutionary code CLES (Code Liegeois d' Evolution Stellaire, 
Scuflaire et al. ,2008aJ . We used the OPAL2001 equa tion o f 
state (Rogers & Nayfonov l2002l Caughlan & Fowler [T988b . 
with nuclear reaction rates from Formicola et al. (I2004I I for the 
^*N{p,yy^O cross-section. Convective transport is treated by 
using the classical mixing length theory of convection (Bohm- 
Vitense 1958) . For the chemical composition, we used the so- 
lar mixture from Asplund et al. (120051) . We used OP opacity 
tables (Seaton I2005I I computed for this mixture. These tables 
are completed at log T < 4. 1 with the low-temperature tables of 
Ferguson et al. (I2005I I. 

When X e [0.682, 0.748 ] and Z e [0.013,0.015] were 
adopted (Przybilla et al. I2008I I. a fine grid of stellar evolution- 
ary models was computed for a mass between 16 and 28 in 
steps of 0. 1 Mq and a core convective overshooting parameter 
Q-Qv between and 0.5 pressure scale heights in steps of 0.05. 
For each main-sequence stellar model, we calculated the theo- 
retical frequency spectrum of low-order p- and g-modes with a 
degree of the oscillation up to ^ = 4 using a standard adiabatic 
code for non-rotating stellar models (Scuflaire et al. i2008bJ . 



5.2. Matching the highest-amplitude frequencies 

To each model of our computed grid of stellar evolutionary mod- 
els, we assigned a ;^'^-value, which compares the observed fre- 
quencies f"^'^ to the theoretically computed ones /J*, as follows: 

/ robs /*th\2 

where cr°''^ is the observational uncertainty. The lower the x^- 
value, the better the match between observed and predicted fre- 
quencies. We selected the observed frequencies with only the 
highest amplitudes, as listed in Table[T] The ;t'^-value is com- 
puted, and the theoretical frequency value is the closest to the 
observed one. Because of the lack of mode identification, no re- 
striction on the degree t was imposed on the theoretical frequen- 
cies, which correspond to zonal modes (m — 0). 

In Fig. Q the positions of all the models having < 2, 
without any a priori restriction on log Teff, but with log g > 
3.89, are compared to the error box representing the position 
of HD 46202 in the log T^ff-logg diagram, as spectroscopically 
deduced. 

We conclude that the positions of our best-fitting models are 
compatible with the spectroscopic one. Indeed, the discrepancy 
found in log g can be considered as small in view of current 
uncertainties in the theory of stellar structure and stellar atmo- 
spheres. Despite great efforts have been made to identify and re- 
duce numerous systematic errors of atmospheric parameters re- 
lated to the data reduction, spectral modelling, and spectral anal- 
ysis (see Nieva & Przvbilla |20 1 Ol for a summary), the remaining 
systematics cannot be excluded. 

One of the main factors in this context is the unknown true 
rotational velocity Viot of HD 46202. The star shows an excep- 
tionally low V sin / for the massive star population of NGC 2244, 
where about 30% of the stars have v sin /> 250 km s"' (Huang 
& Gies2006), indicating a large fraction of truly fast rotators. 
The measured v sin / of HD 46202 is only a lower limit, since 
Viot may reach much higher values if the star is seen nearly pole- 
on. The spectroscopically derived log g would be the polar value 
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Table 3. Physical parameters of HD 46202 as derived from 
the asteroseismic modelling, X e [0.682,0.748] and Z e 
[0.013,0.015]. 



T,s (K) 


34800±780 


log^(cgs) 


3.93±0.03 




24.1±0.8 


Tcvoi (Myr) 


4.3+0.5 


Q'ov 


0.10+0.05 



in that case, while the average logg on the surface of the rota- 
tionally deformed star could be significantly lower. Fast rotation 
would require not only a re-interpretation of the spectroscopi- 
cally derived parameters but also consideration of non-standard 
physics in the stellar structure and evolutionary models (see e.g. 
Ekstrom et al. l2008l) . with additional consequences for the pre- 
dicted frequency spectrum and its excitation. We conclude that 
the potential nature of HD 46202 as a fast rotator seen nearly 
pole-on could resolve the remaining discrepancy between the 
spectroscopic and asteroseismic solutions. The way to check this 
would be to gather a very high-resolution time series (typically 
R~ 100 000 of high S/N (> 300) spectroscopy covering the beat 
pattern of the oscillations). In this way, a spectroscopic mode 
identification can be achieved with certainty and, along with it, 
an estimate of the inclination angle of the star. 

We note that a similar shift in \ogg was also obtained for 
the /? Cep star 6 Ophiuchi (Briquet et al. 12007). For this target, 
the identification of a fundamental radial mode could pinpoint 
the seismic log g-value while the spectroscopic log g-value was 
found to be about 0.15 dex higher. 

5.3. Most probable model parameters 

In Fig.|7] we distinguish two groups of models, shown as open 
squares and as filled circles. Looking deeper into the theoretical 
frequency spectrum of these two groups, we conclude that the 
most probable models are those represented as open squares in 
Fig-E] Our arguments are the following. 

For the models shown as filled circles, the theoretical fre- 
quencies matching the observed frequencies of Table[T] corre- 
spond to /' = 3 or ^ = 4 modes only, except for the frequency /e 
which has an i = 1 . It is very unlikely that no mode with f = 
and ^ = 2 is observed among the highest amplitudes modes. 
The reverse scenario, i.e. with only £ = 0, 1, 2 modes among the 
highest amplitude modes, is certainly more plausible for reasons 
of cancellation effects for higher degree modes. 

The models shown as open squares fit /2 = 4.856 d"' as the 
fundamental radial mode, fj, - 4.691 d"' as an ({,n) = (2, -1) 
mode and /e = 3.004d"' as an {{,n) - (1,-1) mode. For the 
modes with lower frequency values, the theoretical frequency 
spectrum is denser, so an unambiguous identification is more dif- 
ficult. Our best-fitting model has a mass between 23.3 and 24.9 
Mq and a core convective overshooting parameter between 0.05 
and 0.15 pressure scale heights. The other physical parameters 
are given in Table[3] 

5.4. Checking the excitation of the modes 

For all models represented in Fig.|2l we checked the excita- 
tion of the pulsation modes with the linear non-adiabatic code 
MAD developed by Dupret et al. ( I2002I I. We found that none of 



the observed frequencies is theoretically excited using the input 
physics described in Sect. l5.1l 

Such an excitation problem has already been encountered for 
other massive stars. For instance, the two /? Cep stars 12 Lacertae 
(B2III) and v Eridani (B2III) present observed modes that have 
not yet been predicted by pulsation models, although differ- 
ent attempts to solve the discrepancies between theory and ob- 
servations have been suggested in the literature. We refer to 
Dziembowski & Pamyatnykh (120081) . and references therein, for 
the latest discussion of the matter. Another choice of opacity ta- 
ble (e.g. OPAL opacities) and/or metal mixture (e.g. Grevesse 
& Noels 1993 ) does not reconcile the observed and excited 
model frequencies for HD 46202. Miglio et al. ( I2007abl l show 
that more modes are found to be excited in pulsating B-type 
stars when using OP opacities with the Asplund et al. (I2005t 
mixture, as adopted in our study, compared to the case where the 
OPAL opacities are considered. To solve shortcomings in our 
understanding of excitation mechanisms in massive stars and, in 
HD 46202 in particular, is beyond the scope of this paper but de- 
serves to be addressed in depth. One explanation might be that 
current opacities are still underestimated in the region where the 
driving of pulsation modes occurs. 



6. Conclusions 

The space white-light CoRoT photometry reveals the presence 
of /3 Cep-like pulsations of very small amplitudes (~0.1 mmag 
and below) in the 09V star HD 46202. By means of CoRoT 
data, variability associated to stellar pulsations, albeit of a dif- 
ferent nature, was also detected in the O-type binary HD 46149 
(Degroote et al. 120101) and in Plaskett's star (Mahy et al. 2010]). 
These discoveries unambiguously prove there are stellar oscilla- 
tions in O-type stars, while opening the way for promising aster- 
oseismic modelling of this kind of massive star. 

For HD 46202, we found stellar models compatible with 
the observed asteroseismic (pulsation frequencies) and non- 
asteroseismic constraints (effective temperature and surface 
gravity). For the model that best explains the observed oscilla- 
tions, the frequency /2 = 4.856 d ' is identified as the funda- 
mental radial mode. Empirical mode identification is desirable 
in order to confirm this conclusion and to identify the wavenum- 
bers ({, m) of other modes. This task is, however, challenging 
because of the very low amplitude of the modes. 

Another challenge concerns the excitation of the modes. 
Even using OP opacities with the Asplund et al. (120051) solar 
mixture, none of the observed frequencies of HD 46202 is pre- 
dicted to be excited by non-adiabatic computations, which ex- 
plain, at least in general, the excitation of the modes in the less 
massive pulsating B stars. Miglio et al. (2007 ab) extensively ex- 
plored the impact of the choice of opacity tables and metal mix- 
tures on the instability domains for stars with masses between 
2.5 and 18 Mq and for different metallicity values Z. Similar 
studies for higher masses, also including mass loss, are neces- 
sary for a comparison with the newly discovered oscillations in 
massive O-type stars from space data with unprecedented qual- 
ity. Our observational results presented here form a good ba- 
sis for such a theoretical study. Moreover, a detailed exploration 
of the physics in the outer layers leading to the opacity-driven 
modes detected in the 09V star HD 46202 and the modes of 
stochastic nature detected in the 08.5V star HD 46149 will al- 
low an in-depth evaluation of the various excitation mechanisms 
that seem to be at work in slowly rotating O-type stars in the 
main-sequence phase. 
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Table A.l. The frequencies derived from the detrended CoRoT light curve after subsequent prewhitening, their corresponding error, amplitude, 
phase, and signal-to-noise ratio. The signal is written as A, sin[2;r(y;(r - fo) + ^i)]-0.00000169, where to = 2451545 HJD. 



ID 




Error 


filMz) 


A(mag) 


(p(2n rad) 


S/N 


/i 


0.50982 


0.00023 


5.90069 


0.000108 


0.174 


9.28 


/2 


4.85606 


0.00025 


56.20440 


0.000093 


0.087 


11.68 


/3 


4.69107 


0.00029 


54.29479 


0.000080 


-0.292 


10.65 


/4 


2.28961 


0.00029 


26.50012 


0.000071 


0.048 


8.54 


/5 


2.64259 


0.00033 


30.58553 


0.000064 


0.369 


7.69 


/6 


3.00402 


0.00037 


34.76875 


0.000059 


-0.421 


7.23 


fi 


2.19493 


0.00048 


25.40428 


0.000046 


-0.372 


5.72 


h 


1.46182 


0.00051 


16.91921 


0.000043 


0.265 


5.59 


/9 


2.22240 


0.00060 


25.72222 


0.000034 


0.446 


4.84 


/lO 


0.80396 


0.00067 


9.30509 


0.000031 


-0.353 


4.46 


/ll 


2.12660 


0.00071 


24.61343 


0.000031 


0.256 


4.30 


fll 


13.97330 


0.00076 


161.72801 


0.000030 


-0.121 


10.60 


/I3 


2.15488 


0.00078 


24.94074 


0.000028 


-0.072 


3.95 


/l4 


2.35706 


0.00077 


27.28079 


0.000028 


-0.205 


4.09 


/15 


2.49231 


0.00079 


28.84618 


0.000027 


-0.115 


4.09 


/16 


0.54203 


0.00082 


6.27350 


0.000026 


-0.389 


4.00 


fll 


0.57365 


0.00084 


6.63947 


0.000022 


-0.016 


4.03 


/18 


0.32639 


0.00086 


3.77766 


0.000031 


0.392 


3.98 


/l9 


0.40030 


0.00096 


4.63310 


0.000027 


-0.016 


3.61 


fio 


0.22654 


0.00100 


2.62199 


0.000021 


0.148 


3.52 


fll 


0.35462 


0.00101 


4.10440 


0.000023 


0.176 


3.59 


fll 


0.66703 


0.00106 


7.72025 


0.000019 


-0.434 


3.46 


fn 


3.08004 


0.00110 


35.64861 


0.000018 


0.168 


3.39 


fiA 


0.85580 


0.00115 


9.90509 


0.000016 


0.409 


3.31 


fl5 


4.38227 


0.00128 


50.72072 


0.000016 


0.199 


3.02 


fl6 


1.48753 


0.00128 


17.21678 


0.000017 


0.077 


3.06 


fn 


2.00846 


0.00129 


23.24606 


0.000018 


-0.424 


3.07 


fii 


0.99298 


0.00130 


11.49282 


0.000017 


0.275 


3.07 


fl9 


1.11162 


0.00141 


12.86597 


0.000015 


0.072 


2.90 


/3() 


0.71597 


0.00145 


8.28669 


0.000015 


0.496 


2.86 


/31 


2.07063 


0.00145 


23.96562 


0.000015 


0.252 


2.89 


fsi 


0.43954 


0.00147 


5.08727 


0.000015 


-0.155 


2.88 


fi3 


1.67322 


0.00149 


19.36597 


0.000014 


-0.327 


2.88 


/34 


0.25381 


0.00151 


2.93762 


0.000015 


-0.472 


2.90 


/35 


3.02746 


0.00151 


35.04005 


0.000016 


0.275 


2.93 


fi6 


1.91309 


0.00159 


22.14225 


0.000012 


-0.022 


2.83 


fyi 


11.22560 


0.00170 


129.92593 


0.000012 


0.282 


4.65 


f}S 


2.89645 


0.00173 


33.52373 


0.000013 


-0.272 


2.63 


/39 


4.55862 


0.00174 


52.76181 


0.000012 


-0.083 


2.66 


/40 


4.76013 


0.00174 


55.09410 


0.000012 


0.346 


3.10 


/4I 


0.93383 


0.00173 


10.80822 


0.000013 


-0.053 


2.74 


/42 


1.42504 


0.00173 


16.49352 


0.000013 


0.278 


2.77 


/43 


3.38899 


0.00174 


39.22442 


0.000012 


0.381 


2.82 


fu 


0.14792 


0.00177 


1.71204 


0.000012 


-0.493 


2.79 


/45 


2.79809 


0.00179 


32.38530 


0.000011 


0.463 


2.80 


/46 


0.61406 


0.00183 


7.10718 


0.000013 


0.236 


2.78 


/47 


0.78215 


0.00185 


9.05266 


0.000013 


-0.360 


2.79 


/48 


1.76133 


0.00189 


20.38576 


0.000011 


0.025 


2.79 


/49 


1.98728 


0.00191 


23.00093 


0.000013 


0.138 


2.76 


/so 


1.58434 


0.00195 


18.33727 


0.000011 


0.015 


2.77 


/51 


0.95913 


0.00200 


11.10104 


0.000010 


-0.281 


2.73 


/52 


4.43451 


0.00202 


51.32535 


0.000010 


-0.318 


2.83 


/53 


2.32122 


0.00203 


26.86597 


0.000010 


-0.459 


2.78 


/54 


1.13398 


0.00204 


UAlAll 


0.000010 


0.437 


2.80 


/ss 


3.21544 


0.00208 


37.21574 


0.000010 


-0.190 


2.78 




8.77345 


0.00217 


101.54456 


0.000009 


0.458 


3.38 


/57 


1.72569 


0.00219 


19.97326 


0.000010 


-0.045 


2.67 


fss 


1.85355 


0.00213 


21.45313 


0.000010 


-0.492 


2.80 


f59 


0.29858 


0.00216 


3.45579 


0.000010 


0.410 


2.77 
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Table A.l. continued. 



ID 


/(d ') 


Error 




A(mag) 


(piln rad) 


S/N 


/eo 


1.02122 


0.00218 


11.81968 


0.000009 


-0.250 


2.80 


j: 

J61 




U.UUzlO 


1 j.yoooj 


U.UUUUl 1 


U.i /U 


Z.60 


fei 


3.91015 


0.00219 


45.25637 


0.000009 


0.381 


2.74 


/63 


5.91544 


0.00222 


68.46574 


0.000009 


-0.040 


2.89 


/64 


1.24273 


0.00226 


14.38345 


0.000010 


-0.160 


2.85 


/65 


1.95592 


0.00224 


22.63796 


0.000009 


0.479 


2.91 


fee 


1.31797 


0.00228 


15.25428 


0.000009 


-0.195 


2.91 


fei 


0.69020 


0.00223 


7.98843 


0.000009 


0.258 


3.00 



